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Abstract 

Two dimensional square lattice general model of the magnetic dot array is intro- 
duced. In this model the intradot self-energy is predicted via the neural network and 
interdot magnetostatic coupling is approximated by the collection of several dipolar 
terms. The model has been applied to disk-shaped cluster involving 193 ultrathin dots 
and 772 interaction centers. In this case among the intradot magnetic structures re- 
trieved by neural networks the important role play single- vortex magnetization modes. 
Several aspects of the model have been understood numerically by means of the sim- 
ulated annealing method. 

1 Introduction 

In the recent years there is a remarkable progress in the technology of the nanofabrication 
of well defined magnetic materials. The material nanoscience based on the epitaxial and 
lithographic techniques Q allows the fabrication of the regular arrays of the magnetic 
particles-dots of well controlled and interesting shape 0, lattice geometry and composi- 
tion. The increasing technological flexibility calls for further physical ideas, which should 
be incorporated into design of the artificial nanoscale magnetic systems. 

The uniformity of the polarization is the general basic aspect discussed in the con- 
nection with the small magnetic particles. The concept of uniformly polarized particle is 
justified only for the particles of an intermediate size I, |. In the theory [| the magne- 
tostatic coupling was derived for the homogeneously polarized and saturated cylindrical 
dots on a rectangular lattice. More restrictive are conditions of the simulation |^, where 
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each dot of array is substituted by a single dipolar moment. This approximation can be 
used only for monodomain dots separated by a sufficiently large distances. When a dot 
array is represented by a system of the interacting dipoles, the search for the ground state 
configuration leads to the formulation typical for the classical dipolar lattices |0, |8|, ^. 
The violation of the intradot homogeneity stems from the competition between the mag- 
netostatic, anisotropy and exchange energy terms. The analytical model of the dot array 
going towards the non-uniformity was proposed in [10|. In this model the interactions of 
dots were described by the quadrupolar terms. 

The problem of the calculation of the magnetization field of a dot array can be in prin- 
cipal formulated in the terms of classical micromagnetic theory |11]. Due to complexity of 
the problem, the important role in its treatment will play the numerical simulations. They 
require the implementation of the sufficiently dense discretization within the each ferro- 
magnetic dot. As usual, the magnetic part of the system can be subdivided into interacting 



dipoles or grains |12, 13|, small ferromagnetic cubes [0] or finite elements [15[. Then the 
optimum spacing of the mesh nodes is determined by a minimum magnetic length scale 
(exchange, wall) of the system. For the majority of ferromagnetic materials, the com- 
prehensive micromagnetic description is attained when the size of discretization elements 
decreases into the nanometer regime [^]. Thus, the simulation of a single ultrathin dot 
of the micrometer size requires about 10^ nodes, although, the qualitative simulations can 
be realized even for 10^ - 10"^ nodes H, From this we can conclude that both detailed 
and truncated micromagnetic description of many-dot array represents rather demanding 
computational task. We summarize, that principal difficulties of micromagnetic dot array 
analysis come from: (i) the interplay of the phenomena on the intradot (exchange, domain 
wall) and external geometric length scales; (ii) complexity of the magnetic structure of the 
non-uniformly polarized dots; (iii) long-range magnetostatic interdot interactions. 

To make the problem of the magnetic ordering of dot arrays tractable by a moderate 
computer facilities, we have developed method, which works on a much coarser mesh 
than usual discretization schemes allow (except the adaptive and multigrid methods). Its 
general idea is the simultaneous simulation of the intradot - micromagnetic and multidot 
scales. This idea was strongly inspired by the multiscale approach [18|. At the present 
stage of the project, the multidot part of simulation has been developed separately and 
the behaviour of the small-intradot scales has been treated only phenomenologically. The 
approach allows a remarkable increase of the simulation speed, indeed, the price payed 
to the scale separation is the appearance of additional parameters. The completing of 
the project needs support of the algorithms of the parameter estimation developed on the 
basis of the standard micromagnetic simulations. Let us to note that similar problems 
were solved in a cellular automaton version of the molecular dynamics |19]. 

In this paper we presented the results of the simulation of magnetic properties of 
quasi-two-dimensional cluster of ultra-thin magnetic dots on the square lattice. The phe- 
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nomenological aspect of our model is a variable intradot inhomogeneity. The formalism 
we developed for this aim has been adopted from the models of the neural networks. 

Prom the point of view of the information theory, the neural networks are continuous, 
unique mappings constructed from the system of known activation functions. The synaptic 
weights of these activation functions are adjusted by the training process. The standard 
problem, which can be effectively solved by the neural networks is the association of the 
input patterns (in our case inputs are effective magnetic moments) with the desired outputs 
(magnetic self-energy of dot). There are many applications, where neural networks can 
be implemented. They allow interpolation of the data generated by the simulations or 
experiment. The example of the physical application is |2C], where neural network was 



used to fit a complicated analytic potential to the set of ab initio data. In |2l|], the 
Hopfield type of the interaction matrix was suggested to simulate the dynamics of the 
complex protein molecule. The specific magnetic application represents the solution of 
the magnetic inverse problem p^ . 

For the purpose to model the variable magnetic intradot inhomogeneity we adopted the 
theory of the radial basis function networks (RBFN) . The RBFN variant of the neural 
network was chosen, because its ingredient is a straightforward and explicit estimation of 
the synaptic weights, which allows more transparent analysis of the physical symmetries. 

The aim of the paper is to make general presentation of the model and present some 
numerical results. The paper is organized as follows: In Section the model of the dot ar- 



ray energy functional is introduced. This Section consists of two parts: in Subsection 2.1 



we introduced the general phenomenological concept of the nonuniform magnetization. 



which utilizes self-energy interpolation by means of RBFN approach. In Subsection 2.2 
the interaction between dots is introduced. In Section ^ our method is applied to the ul- 
trathin square dots, where tendency of the formation of the vortex intradot phase prevails. 
Section |^ provides some details about the implementation of the simulated annealing al- 
gorithm to the problem of the total energy minimization. Finally, in Section |5| we bring 
examples of the numerical simulations. 

2 Model 

2.1 Intradot self-energy 

The microstate of the system of N magnetic dots is described hy — by — N effective 
magnetic moments 

mi„, i = l,2,...,iV, n = l,2,...,iVc, (1) 
which are associated with the magnetization field of dot M(a;, y) via the volume averages 
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dxdyM(a;,y) . (2) 

(x,y)GAi„ 
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The formal integration is performed here over nth volume element of ith dot labeled 
as Aj„. The saturation magnetic moment mgat = IsVd/^c, where Is is the saturation 
magnetization of the dot, plays role of the normalization factor in Eq.(^. Thus for the 
effective magnetic moments we have the bounding 

mL,x + "^\n,y < 1 > (3) 

where x, y subscripts refer to the Cartesian components of the effective magnetic moments 

where ex and ey are the Cartesian unit vectors. The "softness" 
of niin expressed by Eq.(^ is the important model aspect, which differs from the fun- 
damental Brown's postulate [0]. The reason for this modification is that our model is 
formulated for sufficiently larger elements than classical micromagnetic approach. The 
total magnetization per dot per interaction center is given by 



m 



NN^ . , , 

^ 1 = 1 71=1 



The effective moments characterizing ith dot inside A^-cluster are distributed around the 
dot center Rj and located at Nc positions 

Xj„ = Ri + r„, n = 1,2, . . . ,iVc , (5) 

where r„ are some relative coordinates of the interaction centers. From the assumption 
that identical dots are arranged into array it follows that system of r„ vectors is indepen- 
dent of the dot position inside the cluster. 

In a quasi-two dimensional systems, where mj„ is confined to x — y plane, reduced 
information about ith dot microstate is involved in 2A'c dimensional row vector 

rhi = [ mil , , . . . , irijAr^ ] (6) 



fnil,y,rni2,x,'mi2,y, ■ ■ ■ , rriiM^^x, ''^iN^,y 
center X^i center Xi2 center X^jVc 



(7) 



We continue with the construction of relations associating the effective dot moments 
from Eq.(^ with the corresponding self-energies. Here the intradot self-energy is under- 
stood as a part of the total energy, which includes only anisotropy, exchange and intradot 
magnetostatic energy contributions. The Zeeman term and interdot magnetostatic terms, 
which do not contribute to the self-energy are defined independently of the neural network 
part of the model. 

The construction of the self-energy formula E^'^^^ {•fhi) is based on a proper choice of the 
set of special 2 A'c-dimensional memorized vectors (the input patterns of neural networks) 
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of the type Eq.(|^. The memorized vectors are constructed by putting into row Nc two- 



dimensional vectors pn^ 



Ji) Ai) Jq) Jq) Jq) Ji) 

Pl,x^Pl,y^P2,x^P2,yJ ■ ■ ■ 'PNc,x^PNc,y 



(8) 



Here, the superscript q identifies so cahed feature [^. In our case it is an integer from 
the set Aq = {0, 1, . . . , Q — 1}. In the analogy with Eq.(^), the subscripts oi pn^x,Pnl/ run 
over the interaction centers 

To measure the differences between the configurations we have introduced the Eu- 
clidean norm 

1 1 m 1 1 = rh ■ (9) 

written here for some magnetic moment rh [again encoded via the rule from Eq.(0)]. The 
superscript T from Eq.(^) denotes the vector transposition. 

We start the construction of E'^'^^^ (rhi) by assuming that self-energy is known for Q 
memorized vectors 

i?-lf(p(9)) = , ggAQ, (10) 

where w^"^^ are free parameters of our model. The quality of the interpolation via RBFN 
depends on the choice of the basis functions and corresponding weights. Most convenient 
for our preliminary purposes seems to be the use of Nadaraya- Watson regression estimator 



|24|. According to this, the self-energy input-output relation can be written 



where 



Q-i 

= J2 w^'^^ ip^'^H'm) , (11) 

q=0 

exp f Wp^'^^ - ""T-lP 
= ^ ^— (12) 



are radial basis functions {^'^'^^ {jn)\^ q E Aq} satisfying the normalization conditions 
Q-i 

(rh) = 1 , < ij^'^^m) < 1 . (13) 

The dot index of rh was omitted whenever the distinctions between the individual dots is 
unimportant. The choice of the dispersion d^jax /{V2Q) in Eq.dlD, where 



dmax = max IIp^") -p^"')]! (14) 
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is consistent with the recommendation ||2J]. Speaking in terms of the neural networks 
is the activation function, which determines q—th neuron's response to a given 
input fh and self-energy parameter w^'^\ q € Ag is the optimized weight of the link between 
the input and output layer of the network. 

Having established the formula for the calculation of the self-energy, it is easy to prove 
that if the systems of the memorized vectors is composed from the conjugate vector pairs 
pii) ^ p(i') with the same self-energy parameters 

pig) ^ _p{q') ^ ^(q) ^ ^W) (15) 

RBFN self-energy form including the typical combination of terms 

u;('')(exp(-(Q/(iLx) \\P^'^ " + exp(-(Q/(iLx) I + rh\\')) 

posses the reflection symmetry 

£;^^l^(m) = E'''^^{-m) . (16) 

To analyze the configuration snapshots generated during the simulation process we have 
introduced so called feature map J-. It associates any pattern vector rh with the feature 
q* G Aq I^]. The feature q* identifies the index of a nearest memorized vector p^'^ 

T: q* = T{m), (17) 



11^(9 )_ II = 1 1 p^''^ — m 1 1 . (18) 

The classification of the features performed for the whole cluster gives rise to the A^— component 
vector J^{rhi),J^{rfi2), ■ ■ ■ , J^{mN)- The information about this vector can be concentrated 
to the form of the sample averages 

1 ^ 

i=l 

where 6 is the usual Kronecker symbol. 



2.2 Interdot interactions and interactions with the external field 

The standard assumption about the interdot interactions is that they are essentially mag- 
netostatic [||. By using the concept of the effective moments and interaction centers one 
can construct easily the interdot interaction potential. For this aim we expressed the en- 
ergy contribution E'^^j consisting of the dipole-dipole interactions of the effective moment 
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rrij^ [located at the nth center of ith dot at the position Xj„] with Nc effective moments 
of jth dot (j 7^ i) 

^inj = ^ X! X! ^ ^injs,af3rnjs,l3 , (20) 

a=x,y I3=x,y s=l 



where 



^injs,al3 — Ja/3 I I , -'a/3(.Xj — j— . (Zij 



Because Xj„, X^s are scaled by the lattice spacing a, the energy dimension is absorbed 
into dipolar constant 

A = -^^^^^ , (22) 

where /ir is the relative permeability of the matrix. For the interdot magnetostatic energy 
of dot pair we obtain the expression 

Eff' = E E^nj ■ (23) 
n=l 

The key remaining contribution is the Zeeman energy. For ith dot interacting with the 
external magnetic field H we obtain 

Ef = -msat n -Y^min. (24) 

n=l 

In further, to characterize the external field, we have used the reduced undimensional field 
h = msat H/A, its components hx = h. ■ = hcosO, hy = h ■ ey= hsinO and polar angle 
9. The final form of total energy functional is then given by 

TV f AT 1 

E = J2\ i?-i^(m.) + Ef+ Y: Ef-'' ■ (25) 
i=i [ j=i+i J 

The schematic view on its structure is displayed in Fig.||. According to the scheme of the 
computation, the interacting dots have their counterpart in the interacting RBFN blocks. 
This aspect makes our formulation close to the concept of the interacting neural networks 

n. 
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3 The application to the array of rectangle ultrathin mag- 
netic dots 



In the concrete case we have modelled the ultrathin (quasi-two dimensional) rectangle 
magnetic dots of the square profile x and very small height (Vd/^j) <C Id- Each dot 
is replaced by the four moments A'^c = 4 at the positions r„ = (l/4)£d cos (f (re — 1)) 

+ey sin(|(n — 1)) (see Fig. 2(a)). The magnetic moment of Aj„ element [see Fig.P(b)] 
maintains the saturation value mgat = IsYa/^- After this Eq.(^) modifies to the form 
A = Fd'47(64Wr). 

We have proposed model including Q = 11 memorized patterns (see Table 1) and four 
desired self-energy parameters 

£q for (7 = 0, 

£a for g = 7, 8, 9, 10 . 

The additional parameters of the model are reduction coefficients < 1, < 1, < 1 
(Table 1) introduced to modify the size of memorized effective moments. These co- 
efficients describe the deviations of moments from the saturated value. According to 
Eq.(|2^) one can introduce four subsets of vectors {p*^'^^}, {p^^^p*^^^}, {p'^'^\p'^^\p^^\p^^^^^ 
|p(7r^(8)^p{9)^p(io)|_ "Within to each subset, the vectors correspond to the same self- 
energy. This system includes the vector pairs of the opposite sign p^^^ = —p^^\ p^^^ = 
_p(3)^ p{6) _ — = —p^'^\ p(io) _ _p(8) g^j^^ 2ero memorized vector p^^\ This 
structure guarantee the reflection symmetry of the self-energy given by Eq. (|l^) . 

The exceptional vector p^^^ concerns the integral information from multidomain or 
chaotic magnetization modes of the oscillatory or chaotic character Q. Its occurrence 
is a signature of uncertainty in description of a high momentum magnetization modes. 
Among the patterns memorized and restored by RBFN, we focussed attention to the 
vortex magnetization modes |jl[ 26, 28]. The next two vectors p^^'^ ,p^'^^ encode the 



symmetric vortex and counter- vortex configurations in Fig.y(c), Table 1. Similarly, as in 
the case of p^^^ , the total magnetic moment of the symmetric memorized vortex is zero. For 
the rectangular ultrathin isolated dots of square profile and small crystalline anisotropy. 



the vortex type of magnetic ordering was revealed by the Monte-Carlo simulations |26]. 



This finding was confirmed by the experiments |27, |2^. Vortex modes were also detected 
by the simulations on a cubic particles |jl^ for a weak or zero external magnetic fields. The 
system of vectors {p^^\ P^^\ p^^\ p^^^} belonging to the Stoner-Wohlfart type of single- 
domain particle [|ll| is represented by the four parallel effective moments. This ordering 
can also occur by virtue of the external magnetic or magnetostatic fields. The remaining 
intradot configurations labeled by g = 7, 8, 9, 10 should be the potential sources of (shape) 
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Table 1: The list of eleven memorized configurations of the magnetic moments p^'^^ with 
the corresponding self-energy parameters. The minus sign in p^^) substitutes —1. 



Q 


pi<i) 









• 

• • 
• 


= (0,0,0,0,0,0,0,0) 


So 


1 


= 


1 t 


= K^{0, 1, -,0, 0, -,1, 0) 


£. 


2 




T i 


= Kv(0, -,1, 0, 0, 1, -,0) 




3 


p(3) = 


T 

t t 
1 1 

T 


= Kp (0, 1, 0, 1, 0, 1, 0, 1) 


£p 


4 


= 




= Kp(-,0, -,0, -,0, -,0) 




5 




i 

1 1 
I 


= Kp (0, -,0, -,0, -,0, -) 




6 


= 




= Kp(l, 0, 1, 0, 1, 0, 1, 0) 




7 


= 


T T 


= Kd(0, 1, -,o, 0, 1, -,0) 




8 


p(S) = 


i i 


= KdiO, -,-,0, 0, -,-,0) 




9 


p(9) = 


i i . 


= Kd(0, 0, 0, -,1, 0) 




10 


p(10) = 




= HdiO, 1, 1, 0, 0, 1, 1, 0) 





anisotropy. Let us to note that RBFN approach is not sensitive to the physical nature of 
the anisotropy. 

The principal question arises how to determine seven single-dot parameters k^, k^, 
Kd) "^O) £vi and £d- Further work is needed to combine the present simulations with 
the micromagnetic approaches (see e.g.|l6|, p9| ) incorporating the algorithms of the neu- 
ral network learning [|2^. In this paper the magnetic configurations were selected and 
parametrized in heuristic manner. 



4 The implementation of simulated annealing method 



Simulated annealing jS^] is an optimization technique which operates in a manner anal- 
ogous to the physical process of annealing. In this section we discuss some details of its 
implementation to dot array model. The subject of minimization is energy functional 
Eg .(1251) of the effective magnetic moments. The main parts of adopted algorithm are: 



1. Initial state rhi(t = 0) is generated (or read from the data file). 

2. Cooling schedule. The pseudotemperature T{t) = Tq exp ( —t/to)) relaxes as a func- 
tion of the discrete time t = 0,1, . . . , tmax- 



3.1 



3. The update equation is based on the standard algorithm of Metropolis |31 



Single moment moves. The lattice dot index i £ {1,2, ... , N}, and interaction center 
index nG{l,2,...,A'c} are chosen randomly. At given time t, the current effective 
moment mj„(t) undergoes the stochastic single- moment move m^"^'(t-|-l) = mj„(t)+ 
Vjn^{t) , where Y is two-dimensional vector generated over uniform distribution in 
<— 1,1> X <— 1,1>. The generation is repeated until the constrain |mj„(t) -|- 
?;mY(t)| < 1 is satisfied; is the parameter, which controls the remagnetization 
speed \m{t + 1) — m(t)| < Vm/iN Nc). The encoding of mfjl^^{t) gives rise to the 
row vector m*"^'(t). By taking into account Eq.(25) with the concrete form of the 



Nadaraya- Watson self-energy estimator [Eq.(PT[)] and expression for dipole-dipole 
interaction Eq.(20) we obtained for the energy variation A£'j„, associated with the 
elementary move from mj„(t) to m\^^^{t + 1) as 



Q-i 



AEin = 



w 



q=0 



(27) 



, trial 



a=x,y 



N 

/3=x,y i=ljVi s=l 
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- Acceptance criteria If AEin < 0, the trial configuration is accepted automatically 
and mj„(t + 1) = m*"^'(t + 1). If AEin > 0, the trial configuration is accepted if the 
Boltzmann factor exp(— AE'jn/T) is larger or equal to the random number generated 
uniformly over < 0, 1 >. 

- Quasistatic simulations. In the zero temperature limit, the acceptance criteria re- 
duces to the absolute acceptance if AEin < 0. For this dynamics the energy is mono- 
tonically decreasing function of time. Subsequently, the stochastic motion through 
the phase space tends to the accessible basin of attraction. In the matastable state 
the motion gets stuck for a fixed external magnetic field. The sequence of the 
metastable configurations obtained for a gradually changing external magnetic field 
was used for the calculation of the quasistatic hysteresis loops. 

3.2 Complex intradot moves allows to enhance the effectiveness of annealing and left the 
metastable states. The update according to 3.1 is supplemented by : 

3.2a Interchange moves changing the location of two intradot configurations: rh}^^^'^ = 
fhj, fhj^^^ = rhi belonging to two randomly chosen dots i and j. 

3.2b Feature moves starting with the random choice of the dot index i and its feature 
q G Aq, q / J^^rhi). The suggested move is then given by rhf^^^ = p^^K 

3.2c Reflection moves defined by rhf^^^ = —rhi are of the special importance since 
they conserve the self-energy [see Eq.(|T6|)]. Due to this property the simulated 
system can overcome easily the self-energy barriers. 

For the moves 3.2a, b, c the acceptance probability is given by 



mm < 1, exp — ^ AE'j, 

n=l 



4. Stopping criteria. For t < fmaxj the annealing process follows from the step 2 with 
t <— t + 1 for the move 3.1, or with t <^ t + Nc in the case of the complex intradot 
moves 3.2 a, b, c. 



5 Numerical simulations 

For the model defined in Sect.^ we performed the numerical simulations. We studied finite 
disk-shaped cluster ij^^i^^ter — including N = 193 dots of the size 1^ = a/8. We assume 
that each dot has A'^c = 4 centers at the square lattice [see FigQ. The parameters of 
memorized magnetic configurations are = = 0.95, Kp = 1. Consequently, dmax = 4. 
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Table 2: The numerical tests of Nadaraya- Watson formula calculated for two sets of the 
self-energy parameters. The comparison of desired and retrieved energies [see Eq.(p!o[)]. 
The formula exhibits weakly non-uniform response E^^^^ {p^'^-^'^) to the uniform input 

g{q>l) ^ QOA. 



Q 




^solt(^('?))/A 









20 


27.72 





25.85 


1, 2 





2.84 


60 


56.50 


3, 4, 5, 6 


60 


55.02 


60 


56.88 


7, 8, 9, 10 


40 


40.71 


60 


56.82 



The typical parameters of the simulated annealing have been v^a = 0.2, to = 200iVA'c, 
imax = 5toi and A = Tq < lOA. For the quasistatic simulations we used tmax = lOOA^A'^c- 

The initial simulations were performed for a zero external field and zero self-energy 
parameters. In this case the Monte-Carlo minimization of the energy leads to the rapid fall- 
off of the energy towards the non-colinear antiferromagnetic chains. The configuration is 
displayed in Fig.^(a). From this follows that central part of this cluster corresponds to the 
noncolinear antiferromagnetic phase in agreement with magnetic configuration obtained 
for a system of cylindrical dots and truncated dipolar moments . At the same time the 
surface moments which tend to be parallel to the cluster surface exhibit some kind of the 
frustration |26|. The annealing leads to the ground state estimate E = E^^ = — 38.77AiV. 



The previous value can be understood as the threshold for the competition between 
the interdot and intradot structures. The natural way of the stabilization of the intradot 
vortices is to make the parallel structures of the moments {q = 3,4,5,6) energetically 
unfavorable. In the next we will analyzed more restrictive choice: ?x;('J7^i:2) y \E^^\ with 
the calibration condition t(;('?=i'2) = £^ = Q_ 

The previously simulated system was purely magnetostatic. We follow with the simu- 
lations of the opposite kind of systems, where interdot interactions have been completely 
neglected. For these systems we have constructed a quasi-static hysteresis loops. The re- 
sults have been obtained for the several combinations of the self-energy parameters (com- 
parable with |£'^^|). They are presented in Figs.^(a)-(d). In Table 2 we list differences 
between the desired self-energies and outputs of Nadaraya- Watson estimator. Because 
our parameters are free, the inaccuracy stemming from Nadaraya- Watson formula has no 
principal significance for the quality of the result. In the situations, where precision of 
output becomes to be more relevant, the sofisticated RBFN learning is required [^]. The 
interesting situation has occurred for the self-energy parameters 

£:o = 20A , = , £p = 60A , £d = 40A (28) 

corresponding to Fig.^(a). The field dependence of n(g) observed during the remagnetiza- 
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Table 3: The comparing of the several numerical results obtained the external magnetic 
fields h = 0, 5, 10. Table shows structure of energy contributions (divided by A^A) to the 
ground state. Calculated for the parameters from Eq.(p^). 



energy 


h = 


h = 5 


h = 10 


total 


1.58 


-3.41 


-11.85 


magnetic field 


0.00 


-8.25 


-23.10 


self 


6.37 


13.50 


24.23 


magnetostatic 


-4.78 


-8.67 


-12.98 



tion process for q = 1,2 confirms the vortex stabilization around /i ~ and zero remanence. 
Qualitativelly similar behavior with the vortex annihilation and formation was observed in 
the experimental study |28|. The situation changes dramatically when the strong interdot 
interactions are taken into account. Fig]5| shows that their influence causes the non-zero 
ramanence due to suppression of vortices: n{q) < 0.07 for q = 1,2. More detailed analysis 
of this fact has revealed that surface vortices are more stable than vortices from the central 
zones of the cluster. Two next figures ^(a),(b) show how the anisotropy of the memorized 
configurations is reflected by the hysteresis loops of magnetization components mx{h), 
my{h) constructed for 9 = 20°. 

In the case with the non-zero self-energy parameters interact magnetostatically, the 
choice of the initial conditions of the simulated annealing becomes to be more complicated. 
Several flnal configurations obtained by the annealing process are displayed in Figsj8|(b)- 
(i). The preliminary runs evolving from the initial random state get stuck in the local 
minimum E = 6AXN. The configuration of this metastable state is displayed in Fig.^(b). 
The application of the feature map [see Fig]8|(c)] reveals that the clustering of the dot 
states resembles the formation of the homogeneous domains in Q— state Potts model |33]. 
The lowest energy E = 1.58XN was obtained for the system initialized from the vortex 
state. No essential differences between the pure vortex and mixed vortex-countervortex 
initial conditions were observed, contrary to our expectation evoked by the study of dipolar 
system [^6|. The question of the helicity will require more detailed investigation. Fig.^(d) 
shows that magnetostatic deformation of vortices is rather pronounced feature. After the 
deformation, four-moment vortices acquire the nonzero magnetic moments and resemble 
the fans or vortices with the non-central Neel-type core. Their analysis via J- map shows 
the mixing of vortex features q = 1,2 [ n(l) — n{2) ~ 0.5] and separation of the clusters 
with different vorticity [see FigJ^(d). This aposteriory finding confirms that parameters 
of Eq.(|2^) are sufficient for the stabilization of the vortex ground state for h = 0. Two 
configurations displayed in FigsJ^(e),(f) were obtained for the external magnetic field. 
Their energies are listed in Table 3. For h = hx = 5 the field-deformed vortices resemble 
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the fans ordered into the large-scale wave-like structures. The waves are better visible from 
FigJ^(h) showing the detailed snapshot, where intradot magnetic moments are averaged for 
each dot separatelly. The chaining of field-oriented phase is also visible. For h = = 10 
FigJ^(i) demonstrates the formation of the clusters with q = 6. 

To characterize the anisotropy, we have studied the angular dependence of the magne- 
tization for h = 5 and varying 9 G< 0°, 180° >. In the simulations we have distinguished 
between the clockwise and counter-clockwise field rotation directions. After the anneal- 
ing starting from the purely vortex state mi<i<N{t = 0) = p'^^\ 9 = 0" {9 = 180°) we 
performed the series of quasistatic remagnetization steps at zero temperature. These sim- 
ulations have revealed the hard axes , and easy axes e^; zt By [see Fig.§. In addit ion, 
the model system exhibits the angular hysteresis. These results demonstrate how the out- 
puts of RBFN mimic the biaxial anisotropy and how the anisotropy is reinforced by the 
magnetostatic couplings. 



6 Conclusions 

We believe that very general method we have introduced in this paper will be stimulating 
for the people working in the field of the micromagnetic simulations of the nanoscale 
systems. From the point of view of magnetostatics, direct model improvement is possible 
in many ways: a) near dot interactions can be taken into account more accuratelly by 
including the rectangle-rectangle magnetostatic terms; b) to speed-up the computations 
and to extend the system size one can use the hierarchical summation |13|. We thing that 



more realistic simulations will be possible after the finding of a closer relationship between 
the neural networks and outputs of the standard micromagnetic approaches. 
Acknowledgement: This work was supported by the grant no. 1/6020/99 and by the 
Polish-Slovak international grant. 
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7 Figure Captions 
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Figure 1: The schematic view on the calculation of energy functional using RBFN ap- 
proach: a) the magnetic state of the dot array represented by the system of the eight- 
dimensional vectors; b) RBFN apphcd as a predictor of the self-energy output of ith dot 
from the input vector rhf, c) the interaction of dots (RBFN) mediated by the effective 
magnetic moments. 
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Figure 2: The system of dot centers Ri, R2, . . . , R193 of the simulated cluster. The unit 
vectors Bx and are parallel to the main directions of the square lattice. 
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Figure 3: The basic parameters of the dot array model: (a) ith dot is located at the 
position Rj; the magnetic structure of the dot is represented by the four interaction 
centers and four effective magnetic moments rrij^; (b) two examples of the dots showing 
how the effective magnetic moment arises from the integration of the magnetization field 
over the Aj„ element; (c) the example of the memorized intradot configuration consisting 
of the four magnetic moments. The example of the encoding of intradot configuration by 
means of the eight-dimensional vector p*^^). 
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Figure 4: The hysteresis loops of the isolated dots simulated for RBFN parameters: a) see 
Eq.dH); b) £^ = m\£^ = Q-£^ = m\,£Q = 60A; c) £^ = m\,£^ = 2{)X,£^ = {),£q = 60A; 
d) £^ = 0, £:p = 60A, £a = 60A, £q = Q. 
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Figure 5: The hysteresis loop of the interacting dots in the anisotropic case, calculated for 
9 = 0° and paramters from Eg. ([2^). 
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Figure 7: The angular dependence of the magnetization for interacting (int.) and non- 
interacting (no int.) dot systems confirming the biaxial anisotropy. The arrows indicate 
direction of the field rotation. Calculated for h = 5 and parameters Eq.(28). 
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Figure 8: The system of the low energy configurations revealed by the simulated an- 
nealing process. The case (a) belongs to the zero self-energy parameters and essen- 
tially magnetostatic antiferromagnetic configuration obtained for h = 0. The configu- 
ration shows only the moments obtained for each dot separatelly. The annealing provides 
E = E^^ = -38.77XN. The results obtained in cases (b)-(i) belong to the parame- 
ters from Eg. ([2^). Part (b) displays disordered - quenched metastable microstate of the 
energy E = 6.4AA^ (for h = 0) obtained for the random initial condition and forbiden 
moves 3.2 (a),(b),(c). The metadomain structure was analyzed by J- map [case (c)]. 
The annealing from the initial vortex state stabilizes at the lowest energy E = 1.58XN 
[cases (d),(g)]. Part (g) shows only the dots with q = 1, the remaining part with q = 2 
was removed for the clearness. For the non-zero external fields we constructed snapshots 
(e), (h) = 5, hy = (wave-like ordered structures); (f), (i) hx = 10, hy = ; Part (h) 
shows averaged magnetic moments per each dot. It forms wave-like structure. For (i) the 
formation of ^ = 6 phase clusters is visible (here labels of g = 1, 2 features are removed). 
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